Quantum criticality of a Fermi gas with a spherical dispersion minimum 
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We describe the quantum phase transition of a Fermi gas occurring when the quasiparticle exci- 
tation energy has a minimum in momentum space which crosses zero on a sphere of radius fco 7^ 0. 
The quasiparticles have a universal interaction which controls the physical properties in vicinity 
of the quantum critical point. We discuss possible applications to fermionic superfluids formed by 
pairing two fermion species, near the point where the densities of the two species become unequal. 



The study of quantum and thermal fluctuations near 
second-order, zero temperature quantum phase tran- 
sitions ("quantum criticality") has been an important 
theme in modern condensed matter physics provid- 
ing a unifying framework for many experimental stud- 
ies of heavy fermion compounds and correlated oxides. 
With the advent of experiments on trapped, quantum 
degenerate, ultracold atoms a new arena for the study of 
tunable quantum phase transitions has opened up; the 
most prominent example being the study of superfluid- 
insulator transition of bosonic 87 Rb atoms in an opti- 
cal lattice^. Experiments have also studied the conden- 
sate formed by pairs of fermionic 6 Li atoms with distinct 
hyperfine states, and very recently these have been ex- 
tended to a quantum phase transition into a Fermi liquid 
by unbalancing the density of the hyperfine states [3j . 

This paper will further explore the physics of paired 
fermion systems with unequal densities of the two 
fermion species. Apart from cold atom systems, such 
problems are also of interest to studies of electronic su- 
perconductors in an applied magnetic field and to 
the formation of paired quark condensates at high nu- 
cleon densities |5j. We will introduce and solve a the- 
oretical model of a quantum phase transition that oc- 
curs when the densities of the fermion species is initially 
unbalanced, and describe the universal physical proper- 
ties in its vicinity. The basis of our analysis will be a 
second-order quantum critical point (QCP) which allows 
a systematic renormalization group analysis of arbitrary 
fermion interactions in its vicinity. We will show that 
the QCP is a powerful and unifying vantage point for de- 
scribing a variety of strongly interacting phases obtained 
as the fermion densities and temperature are varied. 

The essential characteristic of our model is that there 
is a conserved U(l) "charge", Q, which is carried by the 
underlying fermions; it is also required that the paired 
fermion condensate is neutral under this U(l) symmetry. 
For the ultracold atom systems, Q is the difference in the 
density of the two fermion species, while for electronic 
superconductors in an applied magnetic field, Q is the 
total spin Sz along the field direction. Now consider the 
situation at zero temperature (T) as a chemical potential 
H conjugate to Q is varied jjj . We are interested in a QCP 
at which (Q) has a non- analytic dependence on /x Q; by a 



(a) 



FIG. 1: Properties of H. (a) For /i < the ground state of 
H is the fermion vacuum, while for /i > the ground state 
has fermions occupying the shaded region with ko — &f < 
k < ko + kp. A is a high momentum cutoff. (6) Crossover 
phase diagram: the dashed lines indicate crossovers between 
the labeled regimes, and the universal interaction applies in 
all three regimes. The full lines are contours of constant (Q). 
In applications to paired fermion problems, all phases are su- 
perfluid, and the FFLO or "breached-pair" phases appear at 
low T in the degenerate fermion regime. Note that for exper- 
iments at fixed (Q) (as in cold atoms) there is a substantial 
intermediate range of T in the quantum critical region, even 
though the ground state is in the degenerate fermion regime. 



shift in /j,, we choose this QCP to be at \i — 0. Typically, 
(Q) is independent of ji for /i < 0, and increases smoothly 
with /i for fi > 0. 

We will describe the situation where for fi < there 
is a gapped fermionic quasiparticle excitation which has 
minimum in its dispersion on a spherical shell of momen- 
tum k with k — ko ^ (the ko = case is considered in 
Ref. 0]). At n = 0, the minimum energy of this fermionic 
excitation vanishes, and the onset in (Q) variation is ac- 
companied by the appearance of two new Fermi surfaces 
at wavevectors below and above ko (see Fig.QJi). One of 
our primary results is that the interactions between low 
lying fermion excitations at wavevectors near k = ko are 
universal, and we will determine some of these interac- 
tions exactly. 

It is worth noting here that solution of the conventional 
weak-coupling Bardeen-Cooper-Schrieffer (BCS) theory 
does indeed yield a minimum in the quasiparticle dis- 
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persion at a fco / 0. However, application to the un- 
balanced fermion case does not yield a QCP as described 
above: instead there is a strong first order transition from 
a state with (Q) = to a state with (Q) ^ at some 
/i < 8J; similar first-order transitions have also been 
obtained in recent mean- field analyses 0, 0, |n|. On 
the other hand, there is recent numerical evidence for 
a continuous transition 12], in a situation with strong 
bare interactions between the fermions and a renormal- 
ized quasiparticle dispersion minimum at a kg ^ 0. We 
will find here that a continuous transition can be present 
for generic short-range interactions between the fermions. 
As is also the case for many other second-order criti- 
cal points [13 > long- wavelength sound- mode fluctuations 
of the compressible environment lead to marginally rel- 
evant flows and a likely weak fluctuation-induced first 
order transition Such a transition is best described 
in the renormalization group framework introduced here, 
and is distinct from that of the weak-coupling mean-field 
analyses. In all cases, the underlying second-order QCP 
remains important for understanding the full scope of the 
crossovers and transitions in the fi-T plane. 

We label the low-lying fermionic excitations in the 
vicinity of k — fco by annihilation operators Ck, so that 
Q = X^k c ik c k; note that the fermions do not carry any 
additional 'flavor' or 'spin' label associated with quantum 
numbers distinct from Q, but it is not difficult to extend 
our analysis to include this. The universal quantum crit- 
ical properties can then be described by the following 
simple Hamiltonian 

ff = 5>(A)-M]cU + 7 £ Vk +q 4- q Ck' Ck , (1) 

k 



k,k',q 



where V q is a generic two-body interaction; A is the vol- 
ume of the system, and the single-particle dispersion is 



;(fc) = (fc-fc ) 2 /(2m*). 



(2) 



For repulsive interactions, H has a QCP at /i = 0, where 
the fermion density increases continuously from zero as /i 
increases. We will determine the critical behavior of the 
system at [i = 0, and also study the properties of the sys- 
tem in the dilute limit kp w y / 2m*|/x| -C fco (see Fig.EJi). 
In two dimensions (2D) we find this low-density phase is a 
crystalline phase with charge density wave (CDW) order; 
this corresponds to the Fulde-Ferrell-Larkin-Ovchinnikov 
(FFLO) stated! G3 

in the fermion pairing problem. On 
the other hand we find in 3D a uniform phase is stable 
against CDW instability when the interaction is suffi- 
ciently weak; this corresponds to the "breached pairing" 
phase [jj. 

The Quantum Critical Point - For simplicity we will 
focus mainly on 2D, and briefly comment on the 3D 
case. Near the QCP at fx = 0, the important low-energy 
fermionic modes are those lie along a shell: fcp — A < fc < 



fco + A, with the cutoff A fc ; these are the modes we 
keep in our model. Thus the phase space structure and 
corresponding kinematic constraints of the present prob- 
lem is quite similar to that of the Fermi liquid fixed point 
studied by Shankar 17] using momentum shell renormal- 
ization group (RG). There are, however, some important 
differences. First, the quadratic fermion dispersion (0) 
changes the scaling dimension of various operators. Sec- 
ond, the ground state at the QCP is known exactly; it is 
simply the fermion vacuum. This immediately leads to 
the exact fermion T — Green's function (which is not 
modified by interactions): 



G(w,k) = l/[iw-e(k)]. 



(3) 



The exact two fermion scattering vertices (or four-point 
functions) r(ki, w\, k' 1; cj^; ki + q, u>\ +to, k^ — q, u>[ — u) 
can be obtained simply by summing ladder diagrams. 
Here we need to divide the scattering vertices into two 
classes, (i) The total momentum of the two particles 
Q| = |ki + k'J 3> A; in this case kinematic constraints 
restrict the scattering processes with momentum transfer 
|q| < A, i.e., only forward scattering is possible in this 
case juj. (ii) |Q| = |ki + kjj <C A; in this case large 
angle scattering is possible, and these are the processes 
in the Cooper channel . 

For case (i) the one-loop correction to T is 



\V(S)f 



•ry , (|ki+q|-fc ) 2 + (|k / 1 -q|-A:o 
2m* 



m*[V{9)] 2 A 2 
" 2tt| sin^ ° g 2m*\n\' 



(4) 



where — uj± + lu[ = lu2 + ^4, 9 is the angle between ki 
and ki, and V{6) = V q=0 - V 2 



2k sin(6>/2) 



The ladder sum 



is sum over a geometric series that yields: 

m*V(6) 



T = F(9M) = V(9)/ 



1 



27r| sin 9 



log 



A 2 



2m*\n\ 



(5) 



Clearly F approaches a universal function in the low- 
energy limit fi — * 0, in that it is independent of the bare 
interaction V , when it is repulsive. 

For case (ii), because k' w — k, it is convenient to 
re-parameterize the vertex as r(ki,wi,— ki + Q,£l — 
uJx;k.2,U2, — k2 + Q, — W2 + O), with Q < A, and de- 
compose it into angular momentum channels: T = 
J2„ U n (Q, n)e in{e ^ B ^) . A similar ladder sum yields 



U n (Q,n)=VJ ll + V n k 



4|ft 



Q 2 



m 

7rA 



where V n 



J-Jd9e- me [V(n -9) 



(6) 

V{6)) (thus n is 



odd), and f(x) is a scaling function with the following 
asymptotic behavior: f[x —* 0) — > 1, and f(x — > oo) — > 

\J log(x). Again the vertex function takes a universal 
form in the limit f2, Q — > for repulsive interactions. 
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It is interesting to interpret the above results in the 
RG langauge. The exact result in Eq. (0 indicates that 
the dynamic exponent z = 2, the correlation length ex- 
ponent v = 1/2, and no fermion anomalous dimension 
rj = 0. However, the presence of the dimensionful scale 
ko means that, strictly speaking, there is no true scale- 
invariant fixed point, and scaling arguments are of lim- 
ited utility. Nevertheless, it is possible to write down RG 
equations and assign scaling dimensions for observables. 
The powers of ko associated with an observable cannot 
be predicted a priori, and explicit calculation is required. 
For forward scattering we obtain the RG equation 



dV{9) 
dlogs 



in ' 



it sin ( 



-[V(9)} 2 , 



(7) 



where s is the momentum re-scaling factor for A. 
This implies forward scattering is marginally irrele- 
vant/relevant for repulsive/ attractive interactions. The 
solution of Eq. (JJJ) also yields the universal structure in 
Eq. In the Cooper channel, we have 



dV n 



to* fen 



V 2 . 



dlogs 7rA(s) ™ 



(8) 



We can remove the explicit dependence of the (5 function 
on A(s) by defining V n = V n /A(s), and the flow equation 
for V n then takes the form 



dV n 
dlog s 



m*k ~ r2 



V 2 . 



(9) 



if it 



We thus find V n flows to a fixed point value of 
is repulsive initially. The physical V n is then obtained by 
rescaling with the observed frequency scale, a conclusion 
consistent with Eq. ©. 

The universal scattering vertices lead to a universal T 
dependence of the quasiparticle scattering rate (or inverse 
lifetime), 1/r, in the quantum-critical regime of Fig. ^p. 
The large density of low energy states associated with the 
quadratic quasiparticle dispersion leads to a short life- 
time, and its value has to be determined by computing 
the scattering rate into quasiparticle states which have 
been self-consistently broadened. A standard computa- 
tion of the scattering cross-section between a quasiparti- 
cle on the Fermi surface and pre-existing thermally ex- 
cited quasiparticle to order T 2 leads to the estimate 



1 

T-2D 



T(m*r) 2 (fc2/ m *)i/2 [Max(Ti i/ T2D )]-y\ (10) 



The only self-consistent solution of this is 

1 C 2 T 2 /3(fc2/ m *)l/3 



T~2D 



Hog A2 I 4 / 3 

I 1U & 2m* T 



(11) 



where C2 is a universal constant of order unity. The scat- 
tering rate is universal in the sense that it is independent 
of the bare interaction. 



If the interaction is attractive, T diverges at low en- 
ergy in all channels, and the strongest divergence is in 
the Cooper channel, signaling the formation of bound 
states before the QCP at p — is reached; such bound 
states appear at p = — V 2 k 2 m* /8, where I is the angular 
momentum channel with the strongest pairing interac- 
tion (i.e., Vi being the most negative). 

Many of these results can be generalized to the 3D 
case with very minor modifications. The most impor- 
tant difference between 2D and 3D is in case (i), the 
forward scattering channel. While in 2D the kinematic 
constraint only allows for small momentum transfer (or 
forward scattering), in 3D the scattering process allows 
for a rotation of the two momenta along the direction 
of Q while preserving the angle between the two mo- 
menta 17], i.e. ki • ki w k 2 • ki or 9^^ « 9^ 2 y 2 « 9. 
The angle of rotation <p = (j)^ — </>k 2 ranges from to 
2tt. Thus in 3D the scattering vertex is parameterized as 
F{9, 0, J7). Performing a Fourier transform with respect 
to 0: F(9,(/),n) = J2n F n( ^)e m4 ', and performing the 
ladder sum yields 

m*k a V n (9) , A 2 



F n (9,n) = V n {9)/ 



1 



4tt|cos(0/2) 



log 



2m*\Vt\ 



(12) 



where 



V n {9) 
with q(9, 1 



d<t>(V q 



V q (8,-K-4,))e 



-in<f> 



(13) 
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2fe o sin(6»/2)sin(0/2). While F„(0,ft) ap- 
proaches a universal function that does not depend on the 
bare interaction in the low energy limit, generically this 
is not the case for F(9, </>, f2). This is because F(9, 0, f2) 
receives contributions from all channels, and the energy 
scale below which F n (9, f2) approaches the universal func- 
tion depends on n; this energy scale approaches zero as 
n — > 00. As a consequence the quasiparticle scattering is 
nonuniversal in 3D; a straightforward calculation similar 
to the 2D case finds that it is linear in T with logarith- 
mic corrections whose form depends on the details of the 
bare interaction. 

The Low-density Phase - We now consider the phase 
with a small [i or fermion density p. In the absence of 
interaction, fermions occupy a thin "Fermi shell" with 
ko — kp < k < fc + kp, where the "Fermi wavevector" 
kp = irp/ko in 2D and w 2 p/kQ in 3D is the half thick- 
ness of the shell (see FigQJi). In the following we show 
that this Fermi shell state has a strong tendency toward 
charge density wave (CDW) ordering, as reflected by en- 
hanced static density susceptibilities at certain wavevec- 
tors; the enhancement is due to the fact that the "Fermi 
velocity" Vp = kp/m* vanishes as kp — > 0, but the size 



of the Fermi surface remains finite and of order k n 



We 



will show that in 2D this enhancement is strong enough 
to render the system unstable against crystallization in 
the low-density limit, in the presence of a weak repulsive 
interactions between the fermions. 
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The strongest enhancement of static density suscepti- 
bility for non-interacting fermions, X (Q), is for Q m 2/cf, 
a nesting vector that connects the inner and outer Fermi 
surfaces that enclose the Fermi shell. A straightforward 
calculation for Q <C ko yields 



o m\ 2w * fc ° y ( Q 
where the dimensionless scaling function 



K 2D (x) 



i r** de 

16n 2 x Jq cos 6 
The result is similar in 3D: 

xUQ) = 2m * k2 ° 



log 



1+x COS ( 



K 3D (x) 



kp 3D \2kp 



1 



1 — x cos t 



Q 



8n 2 x J 





1 


+ xt 


[ dt log 






10 


1 


— xt 



(14) 

(15) 

(16) 
(17) 



Both K2d{x) and Kzd{%) are sharply peaked at x = 1. 
The prefactor k^/kp strongly enhances xo near the QCP; 
this enhancement arises from the nearly perfect nesting 
between Fermi surfaces which are curved in the same 
direction. 

The other nesting wave vectors are Q± = 2/co ± 2kp, 
which connect the opposite ends of the outer and inner 
Fermi surfaces respectively. We hnd 



2m* 
(M 



.1— — / dt 



Mi + f) 

2 + t 2 



(18) 



and X3d(Q±) = m*k /32. We find x°(Q±) diverges in 
the low density limit in 2D while saturates in 3D. 

Now let us go beyond non-interacting fermions, and 
include the interactions computed earlier in a RPA-like 
theory of the susceptibility x- Naively one might expect 
repulsive interactions will immediately lead to a CDW 
instability at the wavevector Q = 2kp due to the di- 
vergent x° as kp — > 0. However, there is a cancelation 
of interaction vertices in this limit for such small mo- 
mentum transfer; as a consequence the strongest CDW 
instability is at larger wavevectors Q± , even though \° 
is less singular there. If we assume the q dependence 
of V q is smooth over the scale of 2k f, which is the case 
for short-range interactions, we can easily sum over the 
particle-hole bubble and ladder diagrams and obtain 



x(Q) 



x°(Q) 



i-{v q=a -v Q )x Q {QY 



(19) 



While x°(Q) i s largest at Q — 2kp, interaction does not 
immediately lead to an instability at this wavevector in 
the low density limit, as for generic short-range interac- 
tions we expect V g —o — VqocQ 2 for small Q. We find in- 
stead that the strongest density instability is at Q = Q±, 
where the interaction effects are stronger. In particular, 



as long as V q=0 - V 2 k > 0, in 2D X2d(Q±) always di- 
verges in the low-density limit kp — ► 0, suggesting that 
the fermions form a charge density wave or crystal state 
in the low density limit. Physically the CDW instability 
at wavevectors that are almost independent of density in 
the low-density limit can be understood in the follow- 
ing manner. In our model the low-energy single fermion 
states are those with k fa k : thus localized states con- 
structed from them must be oscillatory with a character- 
istic wavevector 2ko, similar to the Friedel oscillation; 
this leads to strong CDW instability at such wavevec- 
tors. We further note that in the low-density limit the 
instability has equal strength for Q± = 2kg ± 2kp] as 
a consequence density modulations with characteristic 
wavevectors 0{kp) will be present due to interference. 
We expect their magnitudes to be very large as x is very 
strongly enhanced there for small kp/kg. On the other 
hand the in 3D uniform state appears to be stable against 
weak interaction, due to the fact that X3d(Q±) ^ s finite. 

To summarize, we have presented a new perspective 
on the physics of paired fermion systems, for the case 
in which the fermion densities are unbalanced. It has 
recently become possible to study such ultracold atom 
systems 0, while condensed matter systems have been 
under investigation for some time 0, ■ Essentially all 
previous analyses have been carried out in the framework 
of extensions to BCS theory, seeking to minimize the free 
energy with respect to variations in the magnitude and 
spatial dependence of the order parameter. We have ar- 
gued here that a more rigorous approach is to focus on 
the interactions between the renormalized quasiparticles 
as the system approaches a transition to a gapless phase. 
This allowed us to make a number of controlled predic- 
tions on the influence of strong interactions on the quasi- 
particle dynamics and the structure of proximate phases 
in the fi-T phase diagram. 
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